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Abstract. We report the results from a study of the full sample of ~ 6.031 x 10 7 K — > n^^n 
decays recorded by the NA48/2 experiment at the CERN SPS. As first observed in this ex- 
periment, the invariant mass (Moo) distribution shows a cusp-like anomaly in the region 
around Moo = 2m+, where m+ is the charged pion mass. This anomaly has been interpreted 
as an effect due mainly to the final state charge exchange scattering process tt + tt~ — » 7r°7r° in 
K ± — » 7r ± 7r + 7r _ decay. Fits to the Moo distribution using two different theoretical formulations 
provide the presently most precise determination of ao — <J2, the difference between the tvtv 
S-wave scattering lengths in the isospin I — and 1 — 2 states. Higher-order tvtv rescatter- 
ing terms, included in the two formulations, allow also an independent, though less precise, 
determination of a-2. 



PACS. 13.25.Es Decays of K mesons — 13. 75. Lb Meson-meson interactions — 13.40.Ks Elec- 
tromagnetic corrections to strong- and weak-interaction processes — 14.40.Ao it, K and n 
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Introduction 

The main purpose of the NA48/2 experiment at the CERN 
SPS was to search for direct CP violation in decay to 
three pions PEIS] . The experiment used simultaneous K + 
and K~ beams with momenta of 60 GeV/c propagating 
through the detector along the same beam line. Data were 
collected in 2003-2004, providing large samples of fully re- 
constructed K^- — > 7r ± 7r + 7r~ and if ± — > ^^n decays. 

From the analysis of the data collected in 2003, we have 
already reported the observation of a cusp-like anomaly 
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in the 7r°7r° invariant mass (Moo) distribution of — > 
7r ± 7r°7T decays in the region around Moo = 2m+, where 
m + is the charged pion mass [4 . The existence of this 
threshold anomaly had been first predicted in 1961 by 
Budini and Fonda [5], as a result of the charge exchange 
scattering process n + n~ — > TT°n° in — > 7r ± 7r + 7r _ 
decay. These authors had also suggested that the study 
of this anomaly, once found experimentally, would allow 
the determination of the cross-section for ir + ir~ — » 7r°7r° 
at energies very close to threshold. However, samples of 
— > 7r ± 7r 7r° decay events available in those years were 
not sufficient to observe the effect, nor was the Moo resolu- 
tion. As a consequence, in the absence of any experimental 
verification, the article by Budini and Fonda [5] was for- 
gotten. 

More recently, Cabibbo [5] has proposed an interpreta- 
tion of the cusp-like anomaly along the lines proposed by 
Budini and Fonda [5], but expressing the K ± — ► ^ir ^ 
decay amplitude in terms of the 7r + 7r~ — * 7r°7r° amplitude 
at threshold, a x . In the limit of exact isospin symmetry 
a x can be written as (ao — a2)/3, where ao and ai are 
the S-wave 7T7t scattering lengths in the isospin I — and 
1 = 2 states, respectively. 

Here we report the results from a study of the final 
sample of ~ 6.031 x 10 7 K± -> tt^tt decays. Best fits to 
two independent theoretical formulations of rescattering 
effects in — > n^n ^ and K — > 7r ± 7r + 7r _ decays ([7] 
and [8,9 ) provide a precise determination of ao — a^, and 
an independent, though less precise, determination of a-i. 

1 Beam and detectors 

The layout of the beams and detectors is shown schemati- 
cally in Fig. Q] The two simultaneous beams are produced 
by 400 GeV/c protons impinging on a 40 cm long Be tar- 
get. Particles of opposite charge with a central momen- 
tum of 60 GeV/c and a momentum band of ±3.8% (rms) 
produced at zero angle are selected by two systems of 
dipole magnets forming "achromats" with null total de- 
flection, focusing quadrupoles, muon sweepers and colli- 
mators. With 7x 10 11 protons per pulse of ~ 4.5 s duration 
incident on the target the positive (negative) beam flux at 
the entrance of the decay volume is 3.8 x 10 7 (2.6 x 10 7 ) 
particles per pulse, of which ~ 5.7% (~ 4.9%) are K + 
(K~). The decay volume is a 114 m long vacuum tank 
with a diameter of 1.92 m for the first 66 m, and 2.40 m 
for the rest. 

A detailed description of the detector elements is avail- 
able in [TUJ. Charged particles from decays are mea- 
sured by a magnetic spectrometer consisting of four drift 
chambers (DCH1-DCH4, denoted collectively as DCH) 
and a large-aperture dipole magnet located between DCH2 
and DCH3 [ID] . Each chamber has eight planes of sense 
wires, two horizontal, two vertical and two along each of 
two orthogonal 45° directions. The spectrometer is located 
in a tank filled with helium at atmospheric pressure and 
separated from the decay volume by a thin Kevlar® win- 
dow with a thickess of 0.0031 radiation lengths (Ao). A 
16 cm diameter aluminium vacuum tube centred on the 
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Fig. 1. Schematic side view of the NA48/2 beam line, decay volume and detectors (TAX 17, 18: motorised collimators; 
FDFD/DFDF: focusing quadrupoles; KABES 1-3: beam spectrometer stations (not used in this analysis); DCH1-4: drift cham- 
bers; HOD: scintillator hodoscope; LKr: liquid Krypton calorimeter; HAC: hadron calorimeter; MUV: muon veto). Thick lines 
indicate beam axes, narrow lines indicate the projections of the beam envelopes. Note that the vertical scales are different in 
the left and right part of the figure. 



beam axis runs the length of the spectrometer through 
central holes in the Kevlar window, drift chambers and 
calorimeters. Charged particles are magnetically deflected 
in the horizontal plane by an angle corresponding to a 
transverse momentum kick of 120 MeV/c. The momentum 
resolution of the spectrometer is cr(p)/p = 1.02%© 0.044%p 
(p in GeV/c), as derived from the known properties of 
the spectrometer and checked with the measured invari- 
ant mass resolution of — > tt ± tt + tt~ decays. The mag- 
netic spectrometer is followed by a scintillator hodoscope 
consisting of two planes segmented into horizontal and 
vertical strips and arranged in four quadrants. 

A liquid Krypton calorimeter (LKr) llj is used to re- 
construct 7T° — ► 77 decays. It is an almost homogeneous 
ionization chamber with an active volume of ~ 10 m 3 of 
liquid krypton, segmented transversally into 13248 2 cm 
x 2 cm projective cells by a system of Cu-Be ribbon 
electrodes, and with no longitudinal segmentation. The 
calorimeter is 27 Xq thick and has an energy resolution 
a(E)/E = 0.032/^© 0.09/S 0.0042 (E in GeV). The 
space resolution for single electromagnetic showers can be 
parameterized as a x = a y cm for each 

transverse coordinate x, y. 

An additional hodoscope consisting of a plane of scin- 
tillating fibers is installed in the LKr calorimeter at a 
depth of ~ 9.5 Xq with the purpose of sampling electro- 
magnetic showers. It is divided into four quadrants, each 
consisting of eight bundles of vertical fibers optically con- 
nected to photomultiplier tubes. 



2 Event selection and reconstruction 

The A" ± — * 7r ± 7r°7T decays are selected by a two level 
trigger. The first level requires a signal in at least one 
quadrant of the scintillator hodoscope (Ql) in coincidence 
with the presence of energy depositions in LKr consis- 
tent with at least two photons (NUT). At the second level 
(MBX), an on-line processor receiving the drift chamber 
information reconstructs the momentum of charged parti- 
cles and calculates the missing mass under the assumption 
that the particle is a ^ originating from the decay of a 
60 GeV/c K ± travelling along the nominal beam axis. 
The requirement that the missing mass is not consistent 
with the tt° mass rejects most of the main K — > ^ir 
background. The typical rate of this trigger is ~ 15, 000 
per burst. 

Events with at least one charged particle track having 
a momentum above 5 GeV/c, measured with a maximum 
error of 6% (much larger than the magnetic spectrometer 
resolution), and at least four energy clusters in the LKr, 
each consistent, in terms of size and energy, with the elec- 
tromagnetic shower produced by a photon of energy above 

3 GeV, are selected for further analysis. In addition, the 
relative track and photon timings must be consistent with 
the same event within 10 ns, and the clusters must be in 
time between each other within 5 ns. 

The distance between any two photons in the LKr is 
required to be larger than 10 cm, and the distance between 
each photon and the impact point of any track on the 
LKr front face must exceed 15 cm. Fiducial cuts on the 
distance of each photon from the LKr edges and centre 
are also applied in order to ensure full containment of 
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the electromagnetic showers. In addition, because of the 
presence of ~ 100 LKr cells affected by readout problems 
("dead cells"), the minimum distance between the photon 
and the nearest LKr dead cell is required to be at least 
2 cm. 

At the following step of the analysis we check the con- 
sistency of the surviving events with the K ± — ► ^ir ^ 
decay hypothesis. We assume that each possible pair of 
photons originates from a ir° — ► 77 decay and we calcu- 
late the distance between the ir° decay vertex and the 
LKr front face: 



from the beam axisQ In addition, we require a minimal 
separation between the COG and the charged track coor- 
dinates X t ,Y t , as measured in each DCH plane: 



^COG 



V 2 s nCOG 
1 COG ^max ) 
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COG-track 
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where the limits depend on the COG and track impact 
point distributions at each drift chamber (see Table [T|) . 



Di 
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Table 1. Acceptance cuts on event COG and charged track 
coordinates. 



where E%,Ej are the energies of the i-th and j-th pho- 
ton, respectively, Rij is the distance between their impact 
points on LKr, and mo is the ir° mass. 

Among all possible tt° pairs, only those with values 
differing by less than 500 cm are retained further, and 
the distance D of the decay vertex from the LKr is 
taken as the arithmetic average of the two Dij values. This 
choice gives the best 7r°7r° invariant mass resolution near 
threshold: at Moo = 2to + it is ~ 0.56 MeV/c 2 , increasing 
monotonically to ~ 1.4 MeV/c 2 at the upper edge of the 
physical region. The reconstructed distance of the decay 
vertex from the LKr is further required to be at least 2 m 
downstream of the final beam collimator to exclude 7r°- 
mesons produced from beam particles interacting in the 
collimator material (the downstream end of the final beam 
collimator is at Z = — 18 m). 

Because of the long decay volume, a photon emitted 
at small angle to the beam axis may cross the aluminium 
vacuum tube in the spectrometer or the DCH1 central 
flange, and convert to e + e~ before reaching the LKr. In 
such a case the photon must be rejected because its energy 
cannot be measured precisely. To this purpose, for each 
photon detected in LKr we require that its distance from 
the nominal beam axis at the DCH1 plane must be > 
11 cm, assuming an origin on axis at D — 400 cm. In this 
requirement we take into account the resolution of the D 
measurement (the rms of the difference between D values 
for the two photon pairs distribution is about 180 cm). 

Each surviving ir° pair is then combined with a charged 
particle track, assumed to be a n^. Only those combina- 
tions with a total n^^ir energy between 54 and 66 GeV, 
consistent with the beam energy distribution, are retained, 
and the ir^Tr ^ invariant mass M is calculated, after cor- 
recting the charged track momentum vector for the effect 
of the small measured residual magnetic field in the decay 
volume (this correction uses the decay vertex position, D, 
as obtained from LKr information). 

For each 7r ± 7r°7r° combination, the energy-weighed av- 
erage coordinates (center-of-gravity, COG) XcoGjYcog 
are calculated at each DCH plane using the photon im- 
pact points on LKr and the track parameters measured 
before the magnet (so the event COG is a projection of 
the initial kaon line of flight). Acceptance cuts are then 
applied on the COG radial position on each DCH plane 
in order to select only K ± — > -k^iP-k decays originating 
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The values of i?^° G_tracfc take into account both the 
beam width (the cut is made with respect to each event 
COG rather than to the nominal beam center) and the 
area where the track impact point distribution is still sen- 
sitive to the detailed features of the beam shape. In this 
way the effect of these cuts does not depend strongly on 
the beam shape and on the precise knowledge of the beam 
position in space (during data taking, the average beam 
transverse position was observed to move slightly by up to 
2 mm). This cut removes about 28% of events, mainly at 
large M§ , but the statistical precision of the final results 
on the 7T7T scattering lengths is not affected. 

For events with more than one accepted track-cluster 
combination (~ 1.8% of the total), the K± -> 7r ± 7r°7r° 
decay is selected as the ir^Tr ^ combination minimizing 
a quality estimator based on two variables: the difference 
AD of the two D^ values and the difference AM between 
the 7r ± 7r°7r° invariant mass and the nominal mass |12 : 



/ AD 



( AM 



\rmsD(D) J \rmsM(D) 

where the space and mass resolutions rms d , rms m are 
functions of D, as obtained from the measured AD and 
AM distributions. 

Fig. [21 shows the distribution of AM, the difference 
between the 7r ± 7r°7r° invariant mass and the nominal K 
mass for the selected K ± — > n^ir ^ decays (a total of 
6.031 x 10 7 events). This distribution is dominated by 
the gaussian K ± peak, with a resolution a = 1.3 MeV/c 2 . 
There are small non Gaussian tails originating from uniden- 
tified tt^ — > ^ decay in flight or wrong photon pairing. 
The fraction of events with wrong photon pairing in this 
sample is 0.19%, as estimated by the Monte Carlo simu- 
lation described in the next Section. 



1 The beam is focused at the DCH1 plane, where its width 
s ~ 0.45 cm. 
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Fig. 2. Distribution of the difference between the 7r ± 7r°7r° 
invariant mass and the nominal K mass for the selected 
K — * ^n !: decays. 



Fig. [3]shows the distribution of the square of the 7r°7r° 
invariant mass, Mq , for the final event sample. This distri- 
bution is displayed with a bin width of 0.00015 (GeV/c 2 ) 2 , 
with the 51 s * bin centred at Mq = (2m + ) 2 (for most of 
the physical region the bin width is smaller than the M 2 
resolution, which is 0.00031 (GeV/c 2 ) 2 at M 2 = (2m+) 2 ). 
The cusp at M 2 = (2m+) 2 = 0.07792 (GeV/c 2 ) 2 is clearly 
visible. 



3 Monte Carlo simulation 



Samples of simulated K? 



events ~ 10 times 



larger than the data have been generated using a full de- 
tector simulation based on the GEANT-3 package [13]. 
This Monte Carlo (MC) program takes into account all 
detector effects, including the trigger efficiency and the 
presence of a small number (< 1%) of "dead" LKr cells. 
It also includes the simulation of the beam line; the beam 
parameters are tuned for each SPS burst using fully re- 
constructed — ► 7r ± 7r + 7r _ events, which provide pre- 
cise information on the average beam angles and positions 
with respect to the nominal beam axis. Furthermore, the 
requirement that the average reconstructed tt ± tt + tt^ in- 
variant mass is equal to the nominal K mass for both 
K + and K~ fixes the absolute momentum scale of the 
magnetic spectrometer for each charge sign and magnet 
polarity, and monitors continuously the beam momentum 
distributions during data taking. 

The Dalitz plot distribution of — * tt^^tt decays 
has been generated according to a series expansion in the 
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Fig. 3. a: distribution of Mq , the square of the 7r°7r° invariant 
mass; b: enlargement of a narrow region centred at Mq = 
(2m+) 2 (this point is indicated by the arrow). The statistical 
error bars are also shown in these plots. 



Lorentz-invariant variable u = (S3 — so)/Wi, where Sj = 

(Pk - Pi) 2 (i=l,2,3), s - (si + s 2 + a 3 )/3, Pk {Pi) is the 
K(n) four- momentum, and i = 3 corresponds to the tt^ 



[12] . In our case S3 = Mq Q , and s = (m 2 K 



-2ml- 



4)/3. 

For any given value of the generated 7r°7r° invariant mass 
the simulation provides the detection probability and the 
distribution function for the reconstructed value of Mq . 
This allows the transformation of any theoretical distribu- 
tion into an expected distribution which can be compared 
directly with the measured one. 



4 Determination of the tttt scattering 
lengths a and a 2 

The sudden change of slope ( "cusp" ) observed in the Mg 



distribution at Mq = (2m + ) 2 (see Fig. [3]) can be inter- 
preted [5] [6] as a threshold effect from the decay K ± — > 
7r ± 7r + 7r _ contributing to the K ± — > 7r ± 7r°7r° amplitude 
through the charge exchange reaction 7r + 7r~ — > 7r°7r°. In 
the formulation by Cabibbo [6] the A' ± 
amplitude is described as the sum of two terms 



7T ± 7r 7r decay 



MiK* -> ^±77°^°) = Mo + Mi, 



(1) 



where Mo is the tree level if ± — * 7r ± 7r°7r° weak decay 
amplitude, and Mi is the contribution from the K ± — * 
tt ± tt + 'k~ decay amplitude through ir + n~ — » 7r°7r° charge 
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exchange, with the normalization condition M.\ = at 
Mq = (2m + ) 2 . The contribution M.\ is given by 



M 1 = -2a m m + M + ^(^) -1, (2) 

where a^ is the S-wave 7r + 7r~ charge exchange scatter- 
ing length (threshold amplitude), and M + is the A ± — > 
Tr ± n + n~ decay amplitude at M o = 2m + . .Mi changes 
from real to imaginary at M 00 = 2m + with the conse- 
quence that A^i interferes destructively with A4q in the 
region Moo < 2m + , while it adds quadratically above it. 
In the limit of exact isospin symmetry a x = (ao — fl2)/3, 
where ao and ai are the S-wave 7T7T scattering lengths in 
the 1 = and 1 = 2 states, respectively. 

However, it was shown in ref. [4] that a fit of this simple 
formulation to the NA48/2 Mq distribution in the inter- 
val 0.074 < M 2 < 0.097 (GeV/c 2 ) 2 using a x m+ as a free 
parameter gave only a qualitative description of the data, 
with all data points lying systematically above the fit in 
the region near Mq = (2m + ) 2 . It was also shown in ref. [4] 
that a good fit could be obtained using a more complete 
formulation of 7T7t final state interaction [7] which took 
into account all rescattering processes at the one-loop and 
two-loop level. 

In the following sections we present the determina- 
tion of the 7T7T scattering lengths ao and a,2 by fits of the 
full data set described in Section [2] to two theoretical ap- 
proaches: the Cabibbo-Isidori (CI) formulation [7], and 
the more recent Bern-Bonn (BB) formulation [8|. 

In the CI approach, the structure of the cusp singu- 
larity is treated using unitarity, analiticity and cluster de- 
composition properties of the S'-matrix. The decay ampli- 
tude is expanded in powers of 7T7t scattering lengths up 
to order (scattering length) 2 , and electromagnetic effects 
are omitted. 

The BB approach uses a non-rclativistic Lagrangian 
framework, which automatically satisfies unitarity and 
analiticity constraints, and allows one to include electro- 
magnetic contributions in a standard way [9]. 

In all fits we also need information on the — > 
7r ± 7r + 7r~ decay amplitude. To this purpose, we use a sam- 
ple of 4.709 x 10 8 K ± — > 7r ± 7r + 7r _ decays which are also 
measured in this experiment |14j . 

4.1 Fits using the Cabibbo-Isidori theoretical 
formulation 

In the Cabibbo-Isidori (CI) formulation [7] the weak am- 
plitudes for K ± — > 7r ± 7r°7r° and K ± — > 7r ± 7r + 7r~ decay at 
tree level are written as 

Mo = 1 + ^g Q u + hi u 2 + ifc w 2 , (3) 

M+ = A + {l + ^gu+hu 2 + hv 2 ), (4) 

respectively. In Eq. Q u = (S3 - so)/ m +, where sq = 
(m 2 K + 2toq + m+)/3, while in Eq. (j4j u = (s 3 — s + )/m + , 



where s+ = m 2 K /3 + m 2 ^; for both amplitudes Si = (Pk — 
Pi) 2 , where Pr- (Pi) is the K (ir) four- momentum and i = 
3 corresponds to the odd pion (tv from K — > tt^^tt , 
ir^ from K ± — > 7r ± 7r + 7r _ decay), and v = (s\ — S2) /m 2 + . It 
must be noted that in ref. [7] the v dependence of both am- 
plitudes had been ignored because the coefficients fco and 
fc were consistent with zero from previous experiments. 
Within the very high statistical precision of the present 
experiment this assumption is no longer valid. 

Pion-pion rescattering effects are evaluated by means 
of an expansion in powers of the 7T7t scattering lengths 
around the cusp point, M 2 = (2m + ) 2 . The terms added 
to the tree-level decay matrix elements depend on five S- 
wave scattering lengths which are denoted by a x , a ++ , 
a^ , a+Oi aooi an d describe ir + 7r~ — ► tt° ir° , n + ir + — ► 

7T + 7T + , 7T + 7T~ — > 7T + 7T~ , 7T+7T — > 7T + 7T° , 7T° 7T° — > 7T°7r Scat- 
tering, respectively. In the limit of exact isospin symmetry 
these scattering lengths can all be expressed as linear com- 
binations of ao and 02. 

At tree level, omitting one-photon exchange diagrams, 
isospin symmetry breaking contributions to the elastic nir 
scattering amplitude can be expressed as a function of 
one parameter n = (m 2 + — mfy/m 2 ^ = 0.065 15, 1B1H~7]. In 
particular, the ratio between the threshold amplitudes a x , 
a ++ , a_| , a+Oi aoo and the corresponding isospin sym- 
metric amplitudes - evaluated at the mass - is equal 

tO 1 — n for 7T+7T+ — > 7T+7T+, 7T+7T — > 7T+7T , 7T 7T — > TT°7T° , 

1 + n for 7r + 7r~ — » 7r + 7r~, and 1 + 77/3 for 7r + 7r~ — > 7r°7r°. 
These corrections have been applied in order to extract ao 
and 02 from the fit to the M 2 distribution. 

The CI formulation [7] includes all one-loop and two- 
loop rescattering diagrams and can be used to fit both 
K ± — ► 7r ± 7r°7r° and — > tt ± tt + tt^ decay distributions. 
However, rescattering effects are much smaller in K ± — * 
7r ± 7r + 7r~ than in the if — > 7r ± 7r°7r° decay because the 
invariant mass of any two-pion pair is always > 2m + . In- 
deed, a good fit to the — > tt ± tt + tt^ Dalitz plot [2] can 
be obtained with or without the addition of rescattering 
terms to the tree-level weak amplitude of — > 7r ± 7r + 7r _ 
decay. We have checked that both the values of the best fit 
parameters and their statistical errors, as obtained from 
fits to the Mq distribution of A" ± — > 7r ± 7r°7r° decay, un- 
dergo negligible changes whether or not rescattering ef- 
fects are included in the — > tt^tt^tt - decay amplitude. 
This can be understood from the fact that the K — > 
TT ziz iT + n~ decay amplitude enters into the CI formulation 
of rescattering effects in — > tt^tt !! decays as the 
complete expression given by Eq. ((3]). Thus Eq. with 
parameters extracted from a fit to the — > 7r ± 7r + 7r _ 
data, provides an adequate phenomenological description 
of — > Tr ± TT + TT~ decay which can be used in calculating 
rescattering effects in A" ± — > Tr^TT ^ decay. 

In the fits to the Mqo distribution from A ± — > 7r ± 7r°7r° 
decay, the free parameters are (ao — a2)m+, a2m+, go, ho, 
and an overall normalization constant. The coefficient fco 
cannot be directly obtained from a fit to the Mq distribu- 
tion. Its value is determined independently from the Dalitz 
plot distribution of A ± — > tt^^tt decays, as described 
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in the Appendix. The value fco = 0.0099 is kept fixed in 
the fits. 

All A4+ parameters are fixed from data: the coeffi- 
cients g, h, k are obtained from a separate fit to the 
K ± — > 7r ± 7r + 7r^ decay Dalitz plot [14], using A4+ as 
given by Eq. ((4]), and taking into account Coulomb ef- 
fects; and A + is obtained from the measured ratio, R, 
of the K ± — > tt ± tt + tt~ and K ± — > 7r ± 7r°7r° decay rates, 
R = 3.175 ± 0.050 PH, which is proportional to . The 
fit gives g = -0.2112 ± 0.0002, h = 0.0067 ± 0.0003, k = 
-0.00477 ± 0.00008; and we obtain A + = 1.925 ± 0.015. 
These values are kept fixed in the fits to the Mg distri- 
bution from — > Tr*^ ^ decay. 

As explained in Section[5]all fits are performed over the 
M 2 interval from 0.074094 to 0.104244 (GeV/c 2 ) 2 (bin 
26 to 226). The CI formulation [7] does not include ra- 
diative corrections, which are particularly important near 
Mqq = 2m + , and contribute to the formation of 7r + 7r~ 
atoms ( "pionium" ) . For this reason we first exclude from 
the fit a group of seven consecutive bins centred at AT 2 = 
Am 2 + (an interval of ±0.94 MeV/c 2 in M 00 ). The quality 
of this fit is illustrated in Fig. [4^, which displays the quan- 
tity A = (data - fit)/data as a function of Mq . The small 
excess of events from pionium formation is clearly visible. 



Pionium formation and its dominating decay to 7r°7r° 
are taken into account in the fit by multiplying the content 
of the bin centred at M 2 = 4m+ (bin 51) by 1 + f a tom, 
where 1 + f a tom describes the contribution from pionium 
formation and decay. The pionium width is much narrower 
than the bin width, since its mean lifetime is measured to 
be ~ 3 x 10 -15 s [18]; however, the Mq resolution is taken 
into account in the fits as described in the last paragraph 
of Section 3. The results of a fit with f a tom as a free param- 
eter and with no excluded bins near M§ = 4m^_ are given 
in Tables [3 and [3] (fit CI): the quality of this fit is shown in 
Fig. Hb. The best fit value f atom = 0.0533 ± 0.0091 corre- 
sponds to a rate of K ± — > 7r ± + pionium decay, normalized 
to the K ± -> ti^tt+tt- decay rate, of (1.69 ±0.29) x 10" 5 , 
which is larger than the predicted value ~ 0.8 x 10 -5 [TU 
[213] , As discussed in Section this difference is due to ad- 
ditional radiative effects, which are not taken into account 
in the CI formulation [7] and, contrary to pionium forma- 
tion and decay, affect more than one bin. For this reason 
for the fits without the radiative effects taken into account 
we prefer to fix f a tom = 0.0533 and to exclude from the 
fit the seven consecutive bins centred at Mq = Am 2 ^ . The 
results of this fit are listed as Fit CI a in Tables [2] and [3] 

We have also performed fits using the constraint be- 
tween d2 and a predicted by analyticity and chiral sym- 
metry [5T] (we refer to this constraint as the ChPT con- 
straint): 
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Fig. 4. A = (data - fit) /data versus Mq for the rescattering 
formulation of ref. [7]: a - fit with no pionium formation and 
excluding seven consecutive bins centred at A/q = (2m+) 2 (the 
excluded region is shown by the two vertical dotted lines; b - 
fit with pionium CI (see text). The two vertical dashed lines 
show the Mq interval used in the fit. The point Mq Q = (2m+) 2 
is indicated by the arrow. 



a 2 m+ = (-0.0444 ± 0.0008) -f 
-0.61(a m+ - 0.22) 2 



0.236(a o m + - 0.22) 
-9.9(a o m+-0.22) 3 (5) 



The results of these fits are shown in Tables [2] and 
[3] (fits CI X and CI^). For fit CI X no bins near the cusp 
point are excluded and f a tom is a free parameter, while for 
fit CI\ the seven bins centred at Mq = Am 2 ^ are excluded 
and fatom is kept fixed at the value obtained from fit CI X . 



4.2 Fits using the Bern-Bonn theoretical formulation 

The Bern-Bonn (BB) formulation [8] describes the K — * 
37T decay amplitudes using two expansion parameters: a, 
the generic 7T7t scattering amplitude at threshold; and a 
formal parameter e such that in the K-meson rest frame 
the pion momentum is of order e, and its kinetic energy 
T is of order e 2 . In the formulation of ref. [8] the K — » 
37T decay amplitudes include terms up to 0(e 2 , ae 3 , a 2 e 2 ). 
However, in the formulae used in the fits described below 



a 2 e 2 ). In 



these amplitudes include terms up to 0(e , ae 
the BB formulation the description of the K — > 37r decay 
amplitudes is valid over the full physical region 0. 

At tree level the K — > 37r decay amplitudes are ex- 
pressed as polynomials containing terms in T 3 , T 2 , and 
(Ti — X2) 2 , where T 3 is the kinetic energy of the "odd" 
pion (tt^ from K ± — > ^tt ^, tt^ from K — » 7r ± 7T + 7r _ 
decay) in the rest frame, while T\ and T2 are the 



2 We thank the Bern-Bonn group for providing the computer 
code which calculates the K — > 37r decay amplitudes. 
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kinetic energies of the two same-sign pions. Since these 
variables can be expressed as functions of the relativistic 
invariants u and v defined previously, for consistency with 
the fits described in the previous subsection we prefer to 
use the same forms as given in Eqs. ([3]) and (|3|). It must 
be noted, however, that the best fit polynomial coefficients 
are not expected to be equal to those obtained from the 
fits to the CI formulation 7J because the loop diagram 
contributions are different in the two formulations. 

As for CI, also in the BB formulation rescattering ef- 
fects are much smaller in — > Tr^Tr+vr" than in the 
K ± — ► 7r ± 7T 7r° decay, and a good fit to the M± ± distri- 
bution alone can be obtained with or without the addition 
of rescattering terms to the tree-level weak amplitude of 
— > 7r ± 7r + 7r decay. However, contrary to CI, the coef- 
ficients of the tree-level K — > ■k ± -k + tt~ amplitudes enter 
into the — > 7r ± 7r°7r° rescattering terms in different 
combinations. Therefore, the use of a phenomcnological 
description of the — > 7r ± 7r + 7r~ decay amplitude ex- 
tracted from a fit to if — > 7r ± 7r + 7r~ data alone is not jus- 
tified in this case. Thus, in order to obtain a precision on 
the fit parameters which matches the BB approximation 
level, the value of each coefficient of the — > 7r ± 7r + 7r _ 
tree- level amplitude is obtained from the fit0 

We perform simultaneous fits to two distributions: the 
Mq distribution described in Section [2] and the M±± dis- 
tribution from K ± — > ir ziz Tr + 'K~ decay, obtained as a pro- 
jection of the Dalitz plot described in ref. [14]. This latter 
distribution is made with the same binning as for the M 2 
distribution from K ± — ► Tr ± n n decay and consists of 
4.709 x 10 8 events. 

All fits are performed over the Mq interval from 
0.074094 to 0.104244 (GeV/c 2 ) 2 (bin 26 to 226), and from 
0.080694 to 0.119844 (GeV/c 2 ) 2 (bin 70 to 330) for the 
M± ± distribution from — > 7r ± 7r + 7r _ decay. As for the 
A/qq distribution from K ± — > 7r ± 7r°7r° decay, a very large 
sample of simulated — > tt ± tt + it~ decays (see ref. Q3|) 
is used to obtain the detection probability and the distri- 
bution function for the reconstructed value M±± for any 
generated value of Mf. ± . 

In all fits the free parameters are (ao — a2)ro+ and 
<22m+ (or only aom+ for the fit using the ChPT constraint 
given by Eq. {5J), the coefficients of the tree- level weak 
amplitudes go, ho, g, h, k (see Eqs. (0J[4])), and two over- 
all normalization constants (one for each distribution). 
The coefficient ko (see Eq. (|3J|) is determined indepen- 
dently from a separate fit to the Dalitz plot distribution 
of — > 7r ± 7r°7r° decays (see the Appendix). The fixed 
value fco = 0.0085 is used in the fits. In some of the fits the 
contribution from pionium formation, described by f a tom, 
is also a free parameter. 

Since the detection of — > tt^tt ^ and — > 
7r ± 7r + 7r~ decays involves different detector components 
and different triggers (no use of LKr information is made 

3 Nevertheless, if one fixes the coefficients g, h, k in the fit 
to the values obtained from fits to K — > 7r ± 7r + 7r~ data only 
with or without rescattering terms, the corresponding varia- 
tions of the best fit ao,a2 values are much smaller than the 
ao,ei2 statistical errors. 



to select — > 7r ± 7r + 7r _ decays), the ratio of the de- 
tection efficiencies for the two decay modes is not known 
with the precision needed to extract the value of A + (see 
Eq. from the fit. Therefore, as for the CI fits, also for 
the BB fits A + is obtained from the ratio of the — > 
7r ± 7r + 7r _ and A" ± — >■ n^ir ^ decay rates, measured by 
other experiments, R = 3.175 ± 0.050 [12]. 

Tables [5] and [3] show the results of a fit (fit BB) using 
fatom as a free parameter and including all bins around 
the cusp point in the fit; for fit BBa the value of f a tom is 
fixed and seven bins centred at A/q = 4m + are excluded. 
A comparison with the results of the corresponding CI fits 
(fits CI and CI a, respectively) shows that the difference 
between the best fit values of (a — a 2 )m + is rather small 
(about 3%), while the difference between the two a2m+ 
values is much larger. We note that in the BB fits a2m+ 
has a stronger correlation with other fit parameters than 
in the CI fits (see Tables H and HJ . 

Table 4. Parameter correlations for the CI fits (fit CIa in 
Table 01. 





go ho ao — a^ 02 


go 

h 

ao — ai 

02 


1.000 
-0.701 1.000 

0.777 -0.793 1.000 
-0.902 0.936 -0.869 1.000 



Table 5. Parameter correlations for the BB fits (fit BBa in 
Table 01. 





g<> 


ho 


g 


h 


k ao — ai 


ai 


go 


1.000 












ho 


0.996 


1.000 










g 


-0.970 


-0.960 


1.000 








h 


0.206 


0.181 


-0.247 


1.000 






k 


-0.399 


-0.423 


0.359 


0.803 


1.000 




a — a 2 


-0.853 


-0.817 


0.932 


-0.402 


0.141 1.000 




a~2 


0.976 


0.987 


-0.958 


0.099 


-0.503 -0.794 


1.000 



Fits BB X and BB\ (see Tables and [3J are similar 
to BB and BBa, respectively, but the ChPT constraint 
given by Eq. ((5]) is used. Here the best fit value of aom+ 
agrees well with the value obtained from the CI fit (fit 
CIl). 

5 Radiative effects 

5.1 Radiative correction outside the cusp point 

Radiative corrections to both AT ± — * 7r ± 7r 7r and K ± — * 
7r ± 7r + 7r _ decay channels have been recently studied by ex- 
tending the BB formulation [5j to include real and virtual 
photons [9]. In the if ± rest frame the emission of real 
photons is allowed only for photon energies E < E cut . 

We have performed simultaneous fits to the Mq dis- 
tribution from AT ± — > tt^^tt and to the M±± distribu- 
tion from K^- — > 7r ± 7r + 7r _ decays using the formulation 
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Table 2. Fit results without radiative corrections: nn scattering parameters. Parameter values without errors have been kept 
fixed in the fit or calculated using the constraint between a 2 and ao given by Eq. ([5]). 



Fit 


X Z /NDF 


aom+ 


a,2m + 


(a — a 2 )m+ 


$ atom 


CI 


206.3/195 


0.2334(48) 


-0.0392(80) 


0.2727(46) 


0.0533(91) 


CI A 


201.6/189 


0.2345(50) 


-0.0344(86) 


0.2689(50) 


0.0533 


CI X 


210.6/196 


0.2336(27) 


-0.0413 


0.2749(21) 


0.0441(76) 


cii 


207.6/190 


0.2326(27) 


-0.0415 


0.2741(21) 


0.0441 


BB 


462.9/452 


0.2122(107) 


-0.0693(136) 


0.2815(43) 


0.0530(95) 


BB A 


458.5/446 


0.2182(109) 


-0.0594(143) 


0.2776(48) 


0.0530 


BB X 


467.3/453 


0.2321(33) 


-0.0417 


0.2737(26) 


0.0647(76) 


BB X 


459.8/447 


0.2301(34) 


-0.0421 


0.2722(27) 


0.0647 



Table 3. Fit results without radiative corrections: coefficients of the tree-level K — * 3n weak decay amplitudes. Parameter 
values without errors have been kept fixed in the fit. 



Fit 


go 


ho 


ko 


9 


h 


k 


CI 


0.6512(19) 


-0.0386(23) 


0.0099 


-0.2112 


0.0067 


-0.0048 


CI A 


0.6502(20) 


-0.0375(23) 


0.0099 


-0.2112 


0.0067 


-0.0048 


CI X 


0.6485(9) 


-0.0436(8) 


0.0099 


-0.2112 


0.0067 


-0.0048 


cii 


0.6485(9) 


-0.0438(8) 


0.0099 


-0.2112 


0.0067 


-0.0048 


BB 


0.6117(49) 


-0.0589(56) 


0.0085 


-0.1793(20) 


-0.0015(20) 


-0.0053(23) 


BB A 


0.6154(51) 


-0.0550(57) 


0.0085 


-0.1811(23) 


-0.0012(20) 


-0.0059(22) 


BB X 


0.6215(10) 


-0.0480(9) 


0.0085 


-0.1837(5) 


-0.0011(20) 


-0.0074(20) 


BB X 


0.6215(10) 


-0.0483(9) 


0.0085 


-0.1840(5) 


-0.0008(20) 


-0.0071(20) 



of ref. [9] . Our event selection does not exclude the pres- 
ence of additional photons; however, energetic photons 
emitted in decays result in a reconstructed ^ir ^ 
invariant mass lower than the K mass. We set E cut = 
0.010 GcV in order to be consistent with the measured 
■k^tt^tt invariant mass distribution shown in Fig. [3J (the 
same is true for the ir :iz ir + ir~ invariant mass distribution 
from K ± — ► tt ± tt + tt^ decay measured in this experiment 
[13] ). For each fit we adjust the value of A + (see Eq. ((3])) so 
that the ratio of the -> tt ± tt+tt- and -> 7r ± 7r°7r° 
decay rates is consistent with the measured one [12] . 

The formulation of ref. [§] does not include pionium 
formation, and the — > tt 1 * 1 ^ ^ amplitude, ^00+ j nas 
a non-physical singularity at M§ — (2m + ) 2 . To avoid 
problems in the fits, the square of decay amplitude at the 
center of bin 51, where the singularity occurs, is replaced 
by |A o+| 2 (l + fatom), where A 00+ is the decay amplitude 
of the BB formulation without radiative corrections [5], 
and fatom is again a free parameter. 

The results of simultaneous fits to the M§ distribution 
from K — » 7r ± 7T 7r decays, and to the M±± distribution 
from K ± — > 7r ± 7r + 7r _ decay are shown in Tables [S] and [71 
In all these fits the Mq and M±. ± intervals are equal to 
those of the fits described in Sections 14. Il andl 4. 21 (see Ta- 
bles [5] and [3]) . In fit BB all bins around the cusp point are 
included and f a tom is a free parameter, while in fit BBa 
seven consecutive bins centred at Mq = (2m + ) 2 are ex- 
cluded and fatom is fixed to the value given by fit BB. 
A comparison of fit BB or BBa with radiative correc- 
tions taken into account (Table O with the corresponding 
fits without radiative corrections (fits BB, BBa of Table 
[2]) shows that radiative corrections reduce (ao — a2)m + 



by ~ 9%. However, the change in the best fit value of 
a2?7i + is much larger, possibly suggesting again that the 
determination of this scattering length is affected by large 
theoretical uncertainties. 

Fits BB X and BB\ in Tables M and [7J are similar to 
BB and BBa, respectively, but the constraint between 
ci2 and ao predicted by analyticity and chiral symmetry 
[21"] (see Eq. (JS])) is used. A comparison of fits BB X and 
BB\ with the corresponding fits obtained without radia- 
tive corrections (fits BB X , BB\ of Table shows that 
radiative corrections reduce aom + by ~ 6%. 

For all fits BB X to BB\ in Tables © and [7J the effect of 
changing the maximum allowed photon energy E cut from 
0.005 to 0.020 GeV is found to be negligible. 

No study of radiative corrections has been performed 
in the framework of the CI approach [7J. However, the 
dominating radiative effects (Coulomb interaction and 
photon emission) are independent of the specific approxi- 
mation. Therefore, extracting the relative effect of radia- 
tive corrections from the BB calculation and using it for 
the fit to the CI formula is justified. In order to obtain an 
approximate estimate of radiative effects in this case, we 
have corrected the fit procedure by multiplying the abso- 
lute value of the — > it^it ^ decay amplitude given in 
ref. [7] by IAq^/Aoo+I [55], as obtained in the framework 
of the BB formulation J5[[H] . Because of the non-physical 
singularity of Aq§+ at My — (2m + ) 2 in the BB formula- 
tion, in the calculation of the if — > 7r ± 7r°7r° decay am- 
plitude for the 51 st bin we also multiply the squared am- 
plitude of ref. [7J by 1 + fatom- 

The results of these radiative-corrected fits to the Mq 
distribution from — > 7r ± 7r 7r° decay performed using 
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Table 6. Fit results with electromagnetic corrections: 7T7t scattering parameters. Parameter values without errors have been 
kept fixed in the fit or calculated using the constraint between a 2 and ao given by Eq. 



Fit 


X Z /NDF 




d2m + 


(ao — a 2 )m + 


fatom 


CI 


205.6/195 


0.2391(56) 


-0.0092(91) 


0.2483(45) 


0.0625(92) 


CI A 


202.9/189 


0.2400(59) 


-0.0061(98) 


0.2461(49) 


0.0625 


CI X 


222.1/196 


0.2203(28) 


-0.0443 


0.2646(21) 


0.0420(77) 


cii 


219.7/190 


0.2202(28) 


-0.0444 


0.2645(22) 


0.0420 


BB 


477.4/452 


0.2330(92) 


-0.0241(129) 


0.2571(48) 


0.0631(97) 


BB A 


474.4/446 


0.2350(97) 


-0.0194(140) 


0.2544(53) 


0.0631 


BB X 


479.8/453 


0.2186(32) 


-0.0447 


0.2633(24) 


0.0538(77) 


BB X 


478.1/447 


0.2178(33) 


-0.0449 


0.2627(25) 


0.0538 



Table 7. Fit results with electromagnetic corrections: coefficients of the tree-level K — * 3n weak decay amplitudes. Parameter 
values without errors have been kept fixed in the fit. 



Fit 


5() 


ho 


ko 


9 


h 


k 


CI 


0.6453(22) 


-0.0355(18) 


0.0099 


-0.2112 


0.0067 


-0.0048 


CIa 


0.6446(23) 


-0.0352(18) 


0.0099 


-0.2112 


0.0067 


-0.0048 


CI x 


0.6525(9) 


-0.0433(8) 


0.0099 


-0.2112 


0.0067 


-0.0048 


cii 


0.6526(9) 


-0.0432(8) 


0.0099 


-0.2112 


0.0067 


-0.0048 


BB 


0.6293(47) 


-0.0445(46) 


0.0085 


-0.1928(23) 


-0.0000(20) 


-0.0090(20) 


BB A 


0.6311(51) 


-0.0429(49) 


0.0085 


-0.1938(25) 


0.0004(20) 


-0.0089(20) 


BB X 


0.6219(9) 


-0.0520(9) 


0.0085 


-0.1894(4) 


-0.0003(20) 


-0.0077(19) 


BB X 


0.6220(9) 


-0.0521(9) 


0.0085 


-0.1895(4) 


-0.0002(20) 


-0.0077(19) 



Table 8. Fit parameter correlations for the CI formulation 
with radiative correction (fit CI in Table [B). 





flo 


h 


a - a 2 


0,2 fatom 




1.000 








ho 


-0.629 


1.000 






ao - a 2 


0.794 


-0.719 


1.000 




0,2 


-0.913 


0.883 


-0.873 


1.000 


fatom 


-0.516 


0.387 


-0.650 


0.542 1.000 



Table 9. Fit parameter correlations for the BB formulation 
with radiative correction (fit BB in Table [6}. 





Jo 


ho 


9 


h 


k 




a — a 2 


a 2 


So 


1.000 
















h a 


0.997 


1.000 














g 


-0.972 


-0.965 


1.000 












h 


0.234 


0.220 


-0.255 


1.000 










k 


-0.211 


-0.225 


0.194 


0.889 


1.000 








/„<„,„ 


0.597 


0.570 


-0.652 


0.172 


-0.111 


1.000 






ao — Q2 


-0.870 


-0.843 


0.934 


-0.404 


-0.001 


-0.682 


1.000 




12 


0.977 


0.982 


-0.976 


0.141 


-0.310 


0.597 


-0.839 


1.000 




0.24 



0.25 



0.26 0.27 

(a - o 2 )m+ 



0.28 



0.29 



the CI formula are listed in Tables [6] and [7] (Fits CI to 
CI\). The parameter correlations for two fits which in- 
clude electromagnetic effects are shown in Tables |5] and 

m 

Fig- HI illustrates the fit results for the fits CI and BB 
with and without radiative corrections. All the fits are 
performed using the same — > 7r ± 7r 7r° data sample. 



Fig. 5. 68% confidence level ellipses taking into account the 
statistical uncertainties only. Dashed line ellipses: fits CI and 
BB without radiative corrections. Solid line ellipses: fits CI 
and BB with radiative corrections. The theoretical band al- 
lowed by the ChPT constraint (see Eq. ([SJ)) is shown by the 
dotted curves. 
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5.2 Pionium formation and other electromagnetic 
effects at the cusp point 

Pionium formation in particle decay and in charged parti- 
cle scattering was studied in early theoretical work [201123] , 
but a unified description of its production together with 
other electromagnetic effects near threshold was missing. 

In a more recent approach |24j . electromagnetic ef- 
fects in K ± — > 7T ± 7r 7r decay have been studied in the 
framework of nonrelativistic quantum mechanics using a 
potential model to describe the electromagnetic interac- 
tion between the 7r + 7r _ pair in loop diagrams. This model 
is equivalent to a perturbative one, in which all simple 
sequential 7r + 7r~ loops with electromagnetic interactions 
between the two charged pions are taken into account to 
all orders (including the formation of electromagnetically 
bound final states), but there is no emission of real pho- 
tons and the electromagnetic interaction with the other 
7r ± from the K ± — > 7r ± 7r + 7r" decay is ignored. Because of 
these limitations, the model of ref. [24] cannot be directly 
applied to the full physical region of the K — > tt^tt ^ 
decay; however, contrary to the BB formulation [9 J, its in- 
tegral effect over a narrow region which includes the cusp 
point (Mg = 4m^_) can be calculated. 

We have implemented the electromagnetic effects pre- 
dicted by the model of ref. [24] in the parameterization of 
the CI formulation [7] (the detailed procedure is described 
in Eqs. (6, 7, 8) of ref. [25]). In the theoretical Mq distri- 
bution the electromagnetic correction for the bin centred 
at 4m^_ (bin 51), averaged over the bin, depends on the 
bin width, as it includes contributions from both pionium 
bound states with negligible widths and a very narrow 
peak of unbound ir + ir~ states annihilating to 7r°7r°. For 
the bin width of 0.00015 (GeV/c 2 ) 2 used in the fits, these 
effects increase the content of bin 51 by 5.8%, in agree- 
ment with the results of the fits performed using f a tom as 
a free parameter (see Tables [2] [6} . Thus the model of ref. 
[24 explains why the typical fit result for f a tom is nearly 
twice as large as the prediction for pionium contribution 
only, as calculated in refs. [T9]I20]. 

Near the cusp point the two calculations of electromag- 
netic effects [9] and [241125] are very similar numerically, 
thus increasing the confidence in the central cusp bin ra- 
diative effect calculated using Eq. (8) of ref. [35]. However, 
at larger distances from the cusp the approach of refs. [24] 
[25] leads to deviations from the electromagnetic correc- 
tions of ref. [9] • This can be explained by the fact that the 
model of ref. L 24] takes into account only processes that 
dominate near the cusp point. For this reason we do not 
use this model in the fits, but we consider it as a comple- 
mentary calculation limited to a region very close to the 
cusp point, providing a finite result for the bin centred 
at Mq = Arn 2 ^ which the formulation of ref. [9] does not 
provide. 

6 Systematic uncertainties 

As shown below, all systematic corrections affecting the 
best fit values of the coefficients describing the K ± — > 



^ifi-n® weak amplitude at tree level, go and ho (see Eq. (j^J), 
are found to be much smaller than the statistical errors. 
We use these corrections as additional contributions to the 
systematic uncertainties instead of correcting the central 
values of these parameters. 

For a given fit, we find that the systematic uncertain- 
ties affecting the best fit parameters do not change ap- 
preciably if the fit is performed with or without electro- 
magnetic corrections. In addition, we find that, with the 
exception of f a tom, the systematic uncertainties affecting 
all other parameters are practically the same if in the fit 
the seven consecutive bins centred at Mq = Am 2 + are in- 
cluded (and fatom is used as a free parameter), or if they 
are excluded (and the value of f a tom is fixed). 

For these reasons, we give detailed estimates of the 
systematic uncertainties only for fits CI, CI X , BB, BB X 
performed with the decay amplitude corrected for electro- 
magnetic effects. 

The parameters g, h, k which describe the K ± — ► 
7r~' = 7r^~7r — weak amplitude at tree level are used as free 
parameters when fitting the data to the BB formulation 
[5][2]- However, they enter into the — > ^n ^ de- 
cay amplitude only through rescattering terms, thus we 
do not consider the best fit values of these parameters 
as a measurement of physically important values. Here 
we do not estimate the systematic uncertainties affect- 
ing them and we discuss the uncertainties associated with 
K ± — > 7r ± 7r + 7r~ decay in Section [7] In the study of the 
systematic uncertainties affecting the K ± — > ^ir ^ de- 
cay parameters we fix the values of the K — > 7r ± 7r + 7r _ 
decay parameters g, h, k in the BB formulation to their 
best fit values shown in Table [7] 

The fit interval for the presentation of the final re- 
sults (bins 26-226 of width 0.00015 (GeV/c 2 ) 2 , with bin 
51 centred at 4m 2 + ) has been chosen to minimize the total 
experimental error of the measured ao — a-2- If the upper 
limit of the fit region, s™ ax , is increased, the statistical 
error decreases. All our fits give good x 2 up to rather high 
s max va i ues where the acceptance is small Q However, the 
systematic error increases with s™ ax , especially the contri- 
butions from trigger inefficiency and non-linearity of the 
LKr response. The total experimental error on ao — a,2, 
obtained by adding quadratically the statistical and sys- 
tematic error, has a minimum when the upper limit of the 
fit interval corresponds to bin 226. 



6.1 Acceptance 

The detector acceptance to K ± — > -K^ifiifi decays de- 
pends strongly on the position of the K ± decay vertex 
along the nominal beam axis, Z, so the Z distribution 

4 At the maximum kinematically allowed S3 value the tt^ is 
at rest in the K decay frame. In this case, it moves along the 
K flight path inside the beam vacuum tube and cannot be 
detected. Near this maximum S3 value the acceptance is very 
sensitive to the precise beam shape and position due to the n ± 
narrow angular distribution, and it is difficult to reproduce it 
in the Monte-Carlo simulation. 
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provides a sensitive tool to control the quality of the ac- 
ceptance simulation. 

Fig. [6] shows the comparison between the data and 
Monte-Carlo simulated Z distributions. The small differ- 
ence between the shapes of the two distributions in the 
region Z < disappears when the trigger efficiency cor- 
rection is applied, so this difference is taken into account 
in the contribution to the systematic uncertainties from 
the trigger efficiency (see Tables [TTHT4]) . 

A small difference between the shapes of the two dis- 
tributions is also present in the large Z region in the 
area where the acceptance drops because of the increasing 
probability for the charged pion track to cross the spec- 
trometer too close to the event COG. The effect of this 
acceptance difference has been checked by introducing a 
small mismatch in the track radius cuts between real and 
simulated data, and also by applying small changes to the 
LKr energy scale (equivalent to shifts of the event Z posi- 
tion similar to the effect observed in the acceptance) . The 
corresponding small changes of the fit results are consid- 
ered as the acceptance related contribution to the sys- 
tematic uncertainties (quoted as Acceptance(Z) in Tables 

EHUD. 




-2000 2000 4000 6000 8000 



Z(cm) 

Fig. 6. K — > -k^-k^tt decay Z distributions for data and 
Monte-Carlo simulation, a: Experimental (solid circles) and 
simulated (histogram) distributions, normalized to experimen- 
tal statistics, b: Ratio between the experimental and simu- 
lated distributions. The nominal position of LKr front face is 
at Z = 12108.2 cm. 

The Monte Carlo sample from which the acceptance 
and resolution effects used in the fits are derived, is gen- 
erated under the assumption that the A' ± — * ir^ir ^ 
matrix element, A4, depends only on u. We have stud- 



ied the sensitivity of the fit results to the presence of a 
^-dependent term by adding to lA^ 2 a term of the form 
kov 2 or k' Re(A4)v 2 , consistent with the observed v depen- 
dence in the data. The largest variations of the fit results 
are shown in Tables [TTHT41 as the contributions to the sys- 
tematic uncertainties arising from the simplified matrix 
element used in the Monte Carlo (they are quoted as Ac- 
ceptance (V)). 



6.2 Trigger efficiency 

During data taking in 2003 and 2004 some changes to 
the trigger conditions were introduced following improve- 
ments in detector and electronics performance. In addi- 
tion, different minimum bias triggers with different down- 
scaling factors were used. As a consequence, trigger effects 
have been studied separately for the data samples taken 
during seven periods of uniform trigger conditions. De- 
tails of the trigger efficiency for the A ± — > n^^n decay 
events are given in [TJ[3] . 

As described in Section^ K — > tt^tt ^ events were 
recorded by a first level trigger using signals from the scin- 
tillator hodoscope (Ql) and LKr (NUT), followed by a sec- 
ond level trigger using drift chamber information (MBX) . 
Events were also recorded using other triggers with differ- 
ent downscaling factors for different periods: a minimum 
bias NUT trigger (ignoring both Ql and MBX); and a 
minimum bias Q1*MBX trigger (ignoring LKr informa- 
tion). Using the event samples recorded with these down- 
scaled triggers, and selecting A ± — > n ± Tr t) Tr t) decays as 
described in section [21 it was possible to measure sepa- 
rately two efficiencies: 

1. the efficiency of the minimum bias Q1*MBX trigger 
using the event sample recorded by the minimum bias 
NUT trigger; 

2. the efficiency of the minimum bias NUT trigger using 
the events recorded by the minimum bias Q1*MBX 
trigger. 

These two efficiencies were multiplied together to ob- 
tain the full trigger efficiency. 

The measured efficiencies for seven different periods 
are shown in Fig.[7]as a function of the reconstructed M 2 . 
In the initial data taking periods the samples of minimum 
bias events were rather small, resulting in relatively large 
statistical errors. However, we can improve the estimate 
of the trigger efficiency for these periods under the ad- 
ditional assumption that it is a smooth function of Mq 
(this assumption is justified by the fact that no anomaly 
is expected nor observed in its behaviour). We find that a 
2-nd degree polynomial 

Po + Pi * (M 2 - Am 2 + ) + * (Moo - 4m +) 2 ( 6 ) 

describes well the trigger efficiency over the Mq fit inter- 
val. Moreover, over this interval the dependence is almost 
linear, so we expect a negligible effect on the determina- 
tion of the scattering lengths. 
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Fig. 7. Trigger efficiency e as a function of Mq for the differ- 
ent time periods with different trigger conditions (a-c: 2003, 
d-g: 2004). The errors are defined by the available statistics of 
the event samples recorded by the two minimum bias triggers. 

Table 10. Trigger efficiency corrections for the best fit param- 
eters of fits CI and CI X of Table H 





fit CI 


fit CI X 




0.00056(81) 


0.00111(70) 


h 


0.00136(95) 


0.00136(66) 


(do — a,2)m+ 


-0.00041(67) 




ao?n+ 




0.00065(48) 


a,2m + 


0.00226(190) 




f atom 


0.00070(86) 


-0.00049(82) 



Fits are made separately for each of the data taking 
periods shown in Fig. [71 In a first fit, the Mq Q distribution 
from the data and the corresponding trigger efficiency are 
fitted simultaneously, and the theoretical Mg distribu- 
tion, distorted by the acceptance and resolution effects, is 
multiplied by the corresponding trigger efficiency, as pa- 
rameterized using Eq. ©. The fit to the M§ distribution 
alone is then repeated under the assumption of a fully ef- 
ficient trigger, and the results of the two fits are compared 
to obtain the trigger efficiency correction and its effective 
error. As an example, Table ITD1 lists the trigger corrections 
to the best fit parameters of fits CI and CI X (sec Tabic [6]). 

The trigger corrections are all in agreement with zero 
within their statistical uncertainties. For a conservative 
estimate, we combine in quadrature the corrections and 
their errors to obtain the trigger efficiency contribution 
to the systematic uncertainties of the best fit results (see 
Tables EHU. 



6.3 LKr resolution 

As described in Section [2] the 7r°7r° invariant mass Moo is 
determined using only information from the LKr calorime- 
ter (photon energies and coordinates of their impact points) 
The measurement of the scattering lengths relies, there- 
fore, on the correct description of the Moo resolution in 
the Monte Carlo simulation. 

In order to check the quality of the LKr energy reso- 
lution we cannot use the ir° mass peak in the two-photon 
invariant mass distribution, because the nominal ir° mass 
[12] is used in the reconstruction of the two-photon decay 
vertex (see Section [2]) . We find that a convenient variable 
which is sensitive to all random fluctuations of the LKr 
response, and hence to its energy resolution, is the ra- 
tio m,o /m T o , where m^o and m^o are the measured two- 
photon invariant masses for the more and less energetic 
7T°, respectively, in the same K ± — > 7r ± 7r 7r decay. The 
distributions of this ratio for real and simulated events are 
shown in Fig. [8] One can see that the width of the dis- 
tribution for simulated events is slightly larger than that 
of the data: the rms value of the simulated distribution is 
0.0216, while it is 0.0211 for the data. 

In order to check the sensitivity of the fit results to 
a resolution mismatch of this size, we have smeared the 
measured photon energies in the data by adding a ran- 
dom energy with a Gaussian distribution centred at zero 
and with a = 0.06 GeV (see Fig. [8]). Such a change in- 
creases the rms value of the m,o Im^a distribution from 
0.0211 to 0.0224. A fit is then performed for the data sam- 
ple so modified, and the values of the fit parameters are 
compared with those obtained using no energy smearing. 

The artificial smearing of the photon energies described 
above introduces random shifts of the fit parameters within 
their statistical errors. In order to determine these shifts 
more precisely than allowed by the statistics of a single 
fit, we have repeated the fit eleven times using for each 
fit a data sample obtained by smearing the original pho- 
ton energies with a different series of random numbers, 
as described in the previous paragraph. The shifts of the 
fit parameters, averaged over the eleven fits, represent the 
systematic effects, while the errors on those average val- 
ues are the corresponding uncertainties. Conservatively, 
the quadratic sum of the shifts and their errors is quoted 
as "LKr resolution" in Tables [TT1 [LH 



6.4 LKr non-linearity 

In order to study possible non-linearity effects of the LKr 
calorimeter response to low energy photons, we select ir° 
pairs from — > ^tt ^ events using the following cri- 
teria: 

1. both 7T° — > 77 decays must be close to symmetrical 
(0.45 < ^ < 0.55); 

2. the more energetic 7r° (denoted as it®) must fulfil the 
requirement 

22 GeV < £_o < 26 GeV. 
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tainties in Tables ITTHT4I where they are quoted as "LKr 
non-linearity" . 
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Fig. 8. Distributions of the measured ratio m n o /m n o (see 
text) for the data of 2004. a: solid circles - data events; open 
circles - data events with the LKr cluster energies artificially 
smeared as described in the text; histogram - simulated distri- 
bution, normalized to data statistics, b: corresponding ratios 
of data and simulated distributions. 



For the 7r° pairs selected in such way we define the ratio 
of the two-photon invariant masses, r — m^o/m^o, where 

■k\ is the lower energy tt° . Fig. [9] shows the average ratio 
(r) as a function of E^o/2 for both data and simulated 

events (for symmetric 7r° — * 77 decays E^o jl is the photon 
energy). 

Because of the resolution effects discussed in the pre- 
vious subsectiorjf], (r) depends on the lowest pion energy 
even in the case of perfect LKr linearity. However, as 
shown in Fig. [HI for E n o/2 < 9 GeV the values of (r) 
for simulated events are systematically above those of the 
data, providing evidence for the presence of non-linearity 
effects of the LKr response at low energies. 

To study the importance of these effects, we modify all 
simulated events to account for the observed non-linearity 



multiplying each photon energy by the ratio 



{rpata) 

(rMc) 



, where 



(fData) and {fMc) are the average ratios for data and sim- 
ulated events, respectively. As shown in Fig. [9j the values 
of (r) for the sample of simulated events so modified are 
very close to those of the data. The small shifts of the best 
fit parameters obtained using these non-linearity correc- 
tions are taken as contributions to the systematic uncer- 
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Fig. 9. Average r = m^o/m^o versus E n o/2 for n° pairs from 
K — > n ± TY°TY° decays selected as described in the text. Solid 
circles: data; crosses: simulated events; open circles: simulated 
events corrected for non-linearity (see text). The 7T2 energy is 
divided by 2 to compare with the 7 energy for symmetric 7r° 
decays. 



6.5 Hadronic showers in LKR 

The 7T interaction in the LKr may produce multiple en- 
ergy clusters which are located, in general, near the impact 
point of the ir^ track and in some cases may be identified 
as photons. To reject such "fake" photons a cut on the 
distance d between each photon and the impact point of 
any charged particle track at the LKr front face is imple- 
mented in the event selection, as described in Section 
In order to study the effect of these "fake" photons on the 
best fit parameters we have repeated the fits by varying 
the cut on the distance d between 10 and 25 cm in the se- 
lection of both data and simulated — » tt^^tt events. 
The largest deviations from the results obtained with the 
default cut value (d=15 cm) are taken as contributions to 
the systematic uncertainties (see Tables [TT1 [bH) . 



6.6 Other sources 



5 The small resolution mismatch between data and simulated 
events introduces a negligible effect here. 



The Monte Carlo program includes a complete simula- 
tion of the beam magnet system and collimators with the 
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Fig. 10. Distributions of the reconstructed K momentum 
Pk from the data and from Monte-Carlo simulation (2003 
data), a: solid circles - experimental data; dashed line his- 
togram - simulation; solid line histogram - simulation with 
the corrected K spectrum width, b: corresponding ratios of 
data and simulated spectra. 

purpose of reproducing the correlation between the in- 
cident K momenta and trajectories. However, the ab- 
solute beam momentum scale cannot be modelled with 
the required precision, hence we tune the average value to 
the measured ones for each continuous data taking period 
("run") using — ► 7r ± 7r + 7r~ events which are recorded 
during data taking, and also simulated by the Monte Carlo 
program. 

After this adjustment, a residual systematic difference 
still exists between the measured and simulated K ± mo- 
mentum distributions, as shown in Fig. 1101 In order to 
study the sensitivity of the best fit parameters to this dis- 
tribution, we have corrected the width of the simulated 
K ± momentum distribution to reproduce the measured 
distribution (see Fig. [10]) using a method based on the re- 
jection of simulated events. To minimize the random effect 
of this rejection, a fraction of events has also been removed 
from the uncorrected MC sample in such a way that the 
corrected and uncorrected MC samples have a maximum 
overlap of events and the same statistics. The correspond- 
ing changes of the best fit parameters are included in the 
contributions to the systematic uncertainties and quoted 
as ll P K spectrum" in Tables QIHH 

In order to take into account changes of running condi- 
tions during data taking, the number of simulated — > 
ir^iT ^ events for each run should be proportional to the 
corresponding number of events in the data. However, be- 
cause of changes in the trigger efficiency and in acceptance 



related to minor hardware problems, the ratio between the 
number of simulated and real events varies by a few per- 
cent during the whole data taking period. In order to study 
the effect of the small mismatch between the two samples 
on the best fit parameters, we have made them equal run 
by run by a random rejection of selected events. The corre- 
sponding shifts of the best fit parameters are considered as 
a Monte Carlo time dependent systematic error, and are 
listed in Tables [Tl"HT4l where they are quoted as "MC(T)" . 

Table 11. Fit parameter systematic uncertainties in units of 
10 -4 for the CI formulation with electromagnetic corrections 
(fit CI in Table [6]). The factor m+ which should multiply the 
scattering lengths is omitted for simplicity. 



Source 


9o 


ho 


ao 


Q.2 


ao — a2 


fatom 


Acceptance(Z) 


22 


17 


11 


14 


3 


l 


Acceptance(V) 


9 


3 


5 


6 


1 


3 


Trigger efficiency 


10 


17 


22 


30 


8 


11 


LKr resolution 


4 


2 


11 


17 


7 


56 


LKr nonlinearity 


2 


21 


39 


49 


11 


5 


Pk spectrum 


5 


3 


11 


23 


12 


8 


MC(T) 


3 


2 


4 


1 


5 


25 


ko error 


8 


6 


3 


4 


1 


1 


Hadronic showers 


9 


3 


3 


13 


9 


20 


Total systematic 


29 


33 


49 


67 


22 


66 


Statistical 


22 


18 


56 


92 


45 


93 



Table 12. Fit parameter systematic uncertainties in units of 
10~ 4 for the CI formulation with electromagnetic corrections 
and with the ChPT constraint (fit CI X in Table [6)1. The factor 
m+ which should multiply the scattering lengths is omitted for 
simplicity. 



Source 


90 
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a 2 


ao — a 2 


fatom 


Accept ance(Z) 


24 


14 


4 


1 


3 


9 


Acceptance(V) 


8 


4 


2 





2 





Trigger efficiency 


13 


15 


8 


2 


6 


10 


LKr resolution 





2 


2 





1 


46 


LKr nonlinearity 


12 


13 


13 


3 


10 


31 


Pk spectrum 








2 


1 


2 


5 


MC(T) 


2 


2 


6 


1 


4 


24 


ko error 


7 


7 


1 








2 


Hadronic showers 


5 


3 


4 


1 


3 


19 


Total systematic 


33 


26 


18 


4 


14 


65 


Statistical 


9 


8 


28 


6 


21 


77 



7 External uncertainties 

The most important source of external error is the value 
of |j4+|, obtained from the measured ratio of the if ± — * 
7r ± 7r + 7r~ and K — > ^Tr ^ decay rates, R = 3.175 ± 
0.050 [TJ]. This ratio is proportional to |^4+| 2 , so 

8\A+\/\A+\ = Q.5(8R)/R. 
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Table 13. Fit parameter systematic uncertainties in units of 
1CP 4 for the BB formulation with electromagnetic corrections 
(fit BB in Table [6]). The factor m+ which should multiply the 
scattering lengths is omitted for simplicity. 



Source 


9o 


ho 


a 


a 2 


do - 0,2 


fatom 


Accept ance(Z) 


31 


21 


16 


20 


4 





Acceptance(V) 


6 


1 


7 


8 


1 


4 


Trigger efficiency 


26 


22 


29 


39 


10 


13 


LKr resolution 


10 


9 


21 


29 


9 


60 


LKr nonlinearity 


34 


36 


56 


67 


12 


1 


Pk spectrum 


12 


11 


18 


32 


13 


10 


MC(T) 


2 


1 


4 


1 


5 


25 


ko error 


5 


5 


4 


6 


2 


1 


Hadronic showers 


2 


4 


8 


18 


10 


20 


Total systematic 


56 


50 


72 


94 


25 


70 


Statistical 


47 


46 


92 


129 


48 


97 



8 7T7T scattering lengths: final results 

The BB formulation with radiative corrections [9] provides 
presently the most complete description of rescattering 
effects in K — ► 37r decay. For this reason we use the results 
from the fits to this formulation to present our final results 
on the 7T7T scattering lengths: 

(a - a 2 )m + = 0.2571 ± 0.0048(siai.) 

±0.0025(sj/si.)±0.0014(erf.); (7) 

a 2 m+ = -0.024 ±0M3(stat.) 
±0.009(sysi.) ± 0.002(ert.). (8) 

The values of the 7T7t scattering lengths, (drj — a 2 )m + and 
a2m+, are obtained from fit BB of Tabled In addition 
to the statistical, systematic and external errors discussed 
in the previous sections, these values are affected by a 
theoretical uncertainty. We note that, at the level of ap- 
proximation of the BB and CI amplitude expression used 
in the fits, a difference of 0. 0088(3. 4%)is found between 
the values of (oo — a2)m + and of 0.015(62%) for a 2 m+. 
For the sake of comparison with other independent results 
on the 7T7r scattering lengths we take into account these 
differences as theoretical uncertainty. 

From the measurement of the lifetime of pionium by 
the DIRAC experiment at the CERN PS [15] a value of 
|do — ci2\m + = 0.2641qq2o was deduced which agrees, 
within its quoted uncertainty, with our result (it should be 
noted that this measurement provides only a determina- 
tion of |drj — d2|, while our measurement of K — > 7r ± 7r°7r° 
decay is also sensitive to the sign). 

Previous determinations of the 7T7t scattering lengths 
have also relied on the measurement of K ± — > tt+tt - e v e 
(K e 4) decay. Fig. [ITJ compares our results (Eqs. ([7] [SJ) 
with the results from the most recent analysis of a large 
sample of K e 4 decays, also collected by the NA48/2 col- 
laboration [26] . 

If we use the ChPT constraint (see Eq. ([5])), we obtain 
(see fit BB X of Table© 

(d - a 2 )m+ = 0.2633 ± 0.0024(sia<.) ± 

0.0014(sysi.)±0.0019(ea;t.). (9) 

For this fit the theoretical uncertainty affecting the value 
of ao — a 2 is estimated to be ±2% (±0.0053) from a recent 
study of the effect of adding three-loop diagrams to the 
— > ^Tr ^ decay amplitude [27j in the frame of the 
CI formulation .7] (the goals of this study included a more 
precise estimate of the theoretical uncertainties affecting 
the 7T7T scattering lengths). This theoretical uncertainty 
is smaller than that affecting the result of the fit with 
ao — a 2 and a 2 as free parameters, because the theoreti- 
cal uncertainty on a 2 becomes negligible when using the 
ChPT constraint. 

The 68% confidence level ellipse corresponding to the 
result given by Eq. (J9j) is also shown in Fig. [Tl] together 
with a fit to the K e 4 data which uses the same ChPT 
constraint. The ao — a 2 vs a 2 correlation coefficient for this 



Table 14. Fit parameter systematic uncertainties in units of 
10~ 4 for the BB formulation with electromagnetic corrections 
and with the ChPT constraint (fit BB X in Table©. The factor 
ra+ which should multiply the scattering lengths is omitted for 
simplicity. 



Source 


3o 


h 


a 


a-2 


do - a 2 


fatom 


Acceptance(Z) 


24 


14 


4 


1 


3 


9 


Acceptance(V) 


8 


4 


2 


1 


2 





Trigger efficiency 


14 


16 


9 


2 


7 


8 


LKr resolution 





1 


2 


1 


2 


46 


LKr nonlinearity 


12 


13 


13 


3 


10 


31 


Pk spectrum 








2 


1 


2 


5 


MC(T) 


2 


2 


6 


1 


4 


24 


ko error 


7 


7 











2 


Hadronic showers 


5 


3 


4 


1 


3 


17 


Total systematic 


33 


26 


18 


4 


14 


64 


Statistical 


9 


9 


32 


8 


24 


77 



The typical \A + \ uncertainty is, therefore, S\A + \ w 0.015. 

We have checked the shifts of the fit results due to the 
variation of |^4+| within its uncertainty. Each fit is redone 
twice changing the |^4+| value by +(51^4+ 1 and — <5|A + |. One 
half of the variation of the fit parameters corresponding 
to these two fits is listed in Table [T5l and is taken as the 
external contribution to the full parameter uncertainty. 



Table 15. Contributions to the fit parameter uncertainties (in 
units of 10 -4 ) due to the external error <5|^4+|. 



Fit 


go 


ho 


aom+ 


<J2m + 


(a - a 2 )m + 


fatom 


CI 


3 





27 


14 


13 


l 


CI X 


1 


2 


24 


6 


18 


5 


BB 


5 


3 


32 


18 


14 


1 


BB X 





2 


25 


6 


19 


5 
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Fig. 11. 68% confidence level ellipses corresponding to the 
final results of the present paper (small solid line ellipse: fit 
with the ChPT constraint (see Eq. l[5)l); large solid line ellipse: 
fit using ao — i2 and a 2 as independent parameters) , and from 
K e 4 decay [2H] (small dashed line ellipse: fit with the ChPT 
constraint; large dashed line ellipse: fit using ao and a-z as in- 
dependent parameters). Vertical lines: central value from the 
DIRAC experiment 18! (dotted line) and error limits (dashed 
lines) . The 1-sigma theoretical band allowed by the ChPT con- 
straint (see Eq. |(3J) is shown by the dotted curves. 

figure has been calculated taking into account statistical, 
systematic and external covariances. Its value is —0.774, 
while the statistical correlation alone is —0.839 (see Table 

ED. 



Summary and conclusions 

We have studied the 7r°7r° invariant mass distribution mea- 
sured from the final sample of 6.031 x 10 7 — > tt^^tt 
fully reconstructed decays collected by the NA48/2 exper- 
iment at the CERN SPS. As first observed in this exper- 
iment [J], this distribution shows a cusp-like anomaly at 
Mqo — 2m + which is interpreted as an effect due mainly to 
the final state charge- exchange scattering process ir + n~ — > 
7T 7r° in — > 7r ± 7r + 7r _ decay [MS]. 

Good fits to the Mqo distribution have been obtained 
using two different theoretical formulations [7] and [8,9 , 
all including next-to-leading order rescattering terms. We 
use the results of the fit to the formulation which includes 
radiative corrections [9] to determine the difference ao— 02, 
which enters in the leading-order rescattering term, and 
02) which enters in the higher-order rescattering terms, 
where ao and ai are the 1 = and 1 = 2 S-wave irir 



scattering lengths, respectively. These values are given in 
Eqs. (O and JHJ), while Eq. ((9]) gives the result from a fit 
that uses the constraint between a-i and ao predicted by 
analyticity and chiral symmetry [21] (see Eq. J^J). 

As discussed in Section [5] our results agree with the 
values of the 7T7t scattering lengths obtained from the study 
of K e 4 decay [26j , which have errors of comparable magni- 
tude. The value of a — a 2 as quoted in Eqs. and are 
also in agreement with theoretical calculation performed 
in the framework of Chiral Perturbation Theory [281F29] . 
which predict (ao — 02)771+ = 0.265 ± 0.004. 

We finally note a major difference between — > 
7r ± 7r + 7r~ and K — > 7r ± 7r°7r° decays. In the case of — » 
7r~^~ tt^ tt~ decay there is no cusp singularity in the physical 
region because the invariant mass of any pion pair is al- 
ways > 2m + . As a consequence, rescattering effects can be 
reabsorbed in the values of the Dalitz plot parameters g, 
h, k obtained from fits without rescattering, such as those 
discussed in ref. [14] . On the contrary, a correct description 
of the — > 7T ± 7r°7r Dalitz plot is only possible if rescat- 
tering effects are taken into account to the next-to-leading 
order. Furthermore, the values of the parameters go, ho, 
ko which describe the weak K ± — > 7r ± 7r°7r° amplitude at 
tree level depend on the specific theoretical formulation of 
rescattering effects used to fit the data. 

In a forthcoming paper we propose an empirical pa- 
rameterization capable of giving a description of the 
— > 7r ± 7r°7r° Dalitz plot, which does not rely on any 
7T7T rescattering mechanisms, but nevertheless reproduces 
the cusp anomaly at Moo — 2m+. This parameterization 
is useful for computer simulations of — > 7r ± 7r°7r° decay 
requiring a precise description of all Dalitz plot details. 
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Appendix: Measurement of the k parameter 

In order to measure the kg parameter which describes the 
v 2 dependence of the weak amplitude for A" ± — * 7r ± 7r°7r° 
decay at tree level (see Eq.Q), we have performed fits 
to the tt^tt ?!- Dalitz plot. Because of technical complica- 
tions associated with two-dimensional fits, we do not use 
the results of these fits to determine the scattering lengths, 
but focus mainly on the measurement of fco- 

We use two independent methods. In the first method, 
the Dalitz plot is described by two independent variables: 
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Fig. 12. Projections of the K — > 7r ± 7r°7r Dalitz plot onto 
the cos(#) axis (see text). Full circles: data. Dashed (full) line: 
best fit to the CI formulation [7] with fco = (fco = 0.00974). 

Mg and cos(0), where 9 is the angle between the mo- 
mentum vectors of the and one of the two 7r° in the 
rest frame of the 7r° pair (with this choice of variables 
the Dalitz plot has a rectangular physical boundary) . The 
Mq fit interval is identical to the one used for the one- 
dimensional fits described in Sections 14.11 14.21 but the bin 
width is increased from 0.00015 to 0.0003 (GeV/c 2 ) 2 , and 
four consecutive bins around Mq = Am 2 ^ are excluded. 
The cos(0) variable is divided into 21 equal bins from 
-1.05 to 1.05, but only the interval -0.85 < cos(0) < 0.85 
(17 bins) is used in the fits. 

In order to take into account the distortions of the 
theoretical Dalitz plot due to acceptance and resolution 
effects, a four-dimensional matrix (with dimensions 210 x 
21x210x21)is obtained from the Monte Carlo simulation 
described in Section [3l This matrix is used to transform 
the true simulated Dalitz plot into an expected one which 
can be directly compared with the measured Dalitz plot 
at each step of the \ 2 minimization. 

Fits to the CI formulation [7] are performed with a 
fixed value a-i = —0.044. If the fco parameter is kept fixed 
at zero, the fit quality is very poor (x 2 = 4784.4 for 1237 
degrees of freedom); however, if fc is used as a free param- 
eter in the fit, the best fit value is fc = 0.00974 ±0.00016, 
and x 2 — 1223.5 for 1236 degrees of freedom. The results 
of these two fits are shown in Fig. [12l where the data and 
best fit Dalitz plots are projected onto the cos(0) axis. 

A simultaneous fit to the Dalitz plot from — > 
7r ± 7r°7r° decay and to the M±± distribution from — > 
7r ± 7r + 7r _ decay is performed in the frame of the BB formu- 
lation [5] using the constraint between and ao predicted 
by analyticity and chiral symmetry (see Eq.©). The best 



fit gives fc = 0.00850±0.00014, with x 2 = 1975.5 for 1901 
degrees of freedom. The difference between the fco value so 
obtained and that obtained from a fit to the CI formula- 
tion [7] is due to the rescattering contributions which are 
different in the two formulations. When radiative correc- 
tions are included in the fit [9] , fco is practically unchanged 
(its best fit value is 0.008495), demonstrating that electro- 
magnetic corrections have a negligible effect on its deter- 
mination. 

The second fitting method is based on the event weight- 
ing technique. In order to study the size of the trigger ef- 
fect on the fit parameters, we use a fraction of the data 
taken with uniform trigger conditions and associated with 
a large minimum bias event sample which allows a precise 
evaluation of the trigger efficiency. 

The Dalitz plot is described by the u and \v\ vari- 
ables (see Eq.flHJ), and the intervals —1.45 < u < 1.35 
and \v\ < 2.8 are each sudivided into 50 equal size bins. 
The fits are performed using the CI formulation 7\ over 
a wide region which excludes only the tails of the distri- 
bution (0 < |u| < 0.9 v max , u < 0.9). All bins around 
the cusp point are included, and pionium formation is 
taken into account by multiplying the theoretical — ► 
7r ± 7r°7r° decay probability by the factor 1.055 in the in- 
terval \M§ - 4m 2 _| < 0.000075 (GeV/c 2 ) 2 . The fits are 
performed with a fixed value ai = —0.044. 

In the fits we use the Dalitz plots distributions of the 
selected events, corrected (or not corrected) for the trigger 
efficiency, and of a corresponding subsample of ~ 2.8 x 10 7 
simulated events generated with a simple matrix element 
■M-sim without rescattering effects and with fixed values of 
go, hg and fco- At every iteration in the x 2 minimization, 
each simulated event is reweighted by the ratio r^p^— p, 
where M. is the matrix element which includes rescat- 
tering and is calculated with the new fitting parameters, 
and both M. and M S i m are calculated at the generated 
u, \v\ values. The simulated events so weighted are then 
rebinned, and their two-dimensional u, \v\ distribution is 
compared with that of the data. 

A good fit (x 2 = 1166 for 1257 degrees of freedom) is 
obtained when the trigger efficiency is taken into account, 
giving fc = 0.00966 ± 0.00018. If the trigger effect is ig- 
nored, the x 2 value is somewhat worse (% 2 = 1276) and we 
obtain fc = 0.01010 ± 0.00017. This result demonstrates 
that the trigger effect is important for the wide region of 
the Dalitz plot used in the fit, increasing the measured fco 
by w 0.0004. 

The data used in these fits overlap only partially with 
the data used in the fit to the CI formulation [7] performed 
using the first method and discussed above, but the results 
have almost equal statistical errors. We average the two 
results from the fits without trigger correction, obtaining 
fc = (0.00974 + 0.01010)/2 = 0.0099. We take the sta- 
tistical error of one of them as the statistical error of the 
measured fco value, and conservatively take one half of the 
difference between them as the contribution to the sys- 
tematic error due to the different fitting techniques. As 
mentioned above, the trigger correction shifts the fco cen- 
tral value by —0.0004. Because this effect is measured only 
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with a partial data sample, we also add it in quadrature 
to the systematic error. So our measurement of fco in the 
frame of the CI rescattering formulation [7] gives 

k = 0.0095 ± 0.00017 (stat.) ± 0.00048(syst.) 

= 0.0095 ±0.0005. 

For most of the one-dimensional fits discussed in the 
present paper we do not apply any trigger correction, so 
here we use the effective value ko = 0.0099 for the fits 
to the CI formulation [7], and fco = 0.0085 for the fits to 
the BB formulation [ElE]. Since ko is kept fixed in those 
fits, we check the variations of all the best fit parameters 
by varying fco within the limits defined by its full error. 
These variations are listed in Tables [TTHT41 where they are 
denoted as "fco error" . 
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